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We address the question of the quantitative relationship between thermodynamic phase transitions 
and topological changes in the potential energy manifold analyzing two classes of one dimensional 
models, the Burkhardt solid-on-solid model and the Peyrard-Bishop model for DNA thermal de- 
naturation, both in the confining and non-confining version. These models, apparently, do not fit 
[M. Kastner, Phys. Rev. Lett. 93, 150601 (2004)] in the general idea that the phase transition is 
signaled by a topological discontinuity. We show that in both models the phase transition energy 
v c is actually non-coincident with, and always higher than, the energy vg at which a topological 
change appears. However, applying a procedure already successfully employed in other cases as the 
mean field 4 model, i. e. introducing a map M : v — *• v s from levels of the energy hypersurface V 
to the level of the stationary points "visited" at temperature T, we find that M(v c )—ve. This result 
enhances the relevance of the underlying stationary points in determining the thermodynamics of a 
system, and extends the validity of the topological approach to the study of phase transition to the 
elusive one-dimensional systems considered here. 

PACS numbers: 64.60.-i, 05.70.Fh, 02.40.-k 



Phase transitions are a very well understood subject 
in statistical mechanics. They have been characterized 
in many different ways in the last century and many ob- 
servables related to the phase transition (e.g. critical ex- 
ponents, correlation lengths, etc.) have been computed 
and measured with very high accuracy 0]. 

Recently, a new characterization of phase transitions 
has been proposed by Pettini and coworkers @, 0, 0, 
5]. These authors conjectured that, for classical sys- 
tems defined by a continuous potential energy function 
V{{li}i=i...N), a thermodynamic phase transition, oc- 
curring at a temperature T c , is the manifestation of a 
topological discontinuity, taking place at a specific value 
v(T c ) = N^ 1 (V(q)) (where (•) is the statistical aver- 
age at temperature T) of the potential energy function 
V, or, more precisely, taking place on the hypersurface 
E„ = «y) e R N \ V(qi,...,q N ) = Nv}, at 

v = v(T c ). The most striking consequence of this hypoth- 
esis is that the signature of a phase transition is present 
in the topology of the configuration space independently 
on the statistical measure defined on it. 

The changes in the topology are identified through the 
Morse theory [|| : according to this theory the topological 
changes in a manifold like £„ are related to the presence 
of stationary points of V (points for which VV=0) at 
energy v. However, the precise meaning of the correla- 
tion between topological changes and phase transitions in 
the general case is still a open question. From one hand 
there is a theorem of Franzosi and Pettini _5j , asserting 
that, for "smooth, finite-range and confining microscopic 
interaction potentials V with continuously varying coor- 
dinates,. . . , a topology change of the {S^j^gR at some 
vg is a necessary condition for a phase transition to take 



place at the corresponding energy . . . value" p| . On the 
other hand, there are different numerical studies of vari- 
ous models 

SHSHGUIHEEH (almost all with 
potentials V not fulfilling the hypotheses of the theorem) 
for which a variety of results has been obtained: some are 
in agreement with the "topological hypothesis" of Pettini 
and coworkers, others seem to indicate its failure. 

It is important to underline that the theorem in Ref. [{| 
establishes a necessary condition for a phase transition 
to take place. The problem to find sufficiency condi- 
tions is still an open problem, as pointed out also by 
the proponents of the hypothesis. This question has 
been addressed in two recent papers: in Ref. a one- 
dimensional model, the Burkhardt model with noncon- 
fining potential (see below) was investigated, finding that 
a topological change is present without a phase transition 
at finite temperature; in Ref. [Til ] the mean field spher- 
ical ferromagnet was considered, and it was found that 
the same topological changes happen either in absence 
and in presence of a magnetic field, while in the latter 
case no phase transition occurs (it is worth noting, how- 
ever, that these two models do not fulfill the hypotheses 
of the theorem in Ref. 5]). Moreover, in the thermo- 
dynamic limit, it is likely that, for any finite interval of 
energy / = [ui,V2], there is always a stationary point 
of V(q) with energy v G /; thus, in the thermodynamic 
limit a topological change occurs with probability one in 
any finite interval of energy. Most of these topological 
changes are -obviously- not related to phase transitions, 
and indeed it seems that a topological change must be 
strong enough to be related to a thermodynamic phase 
transition (see Ref. Q for a detailed discussion of this 
point). An indication coming from the analysis of the 
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models cited above is that the presence of a phase tran- 
sition should be related to the presence of a singularity 
in the Euler characteristic at a given energy vg. Thus, 
basically, the idea of the "topological approach is that 
the phase transitions are correlated to abrupt changes in 
topological quantities defined on the stationary points of 
V (as, for example, the Euler characteristic). 

Another open question is the equivalence between the 
energies at which phase transition (v c ) and change in 
the topology (vg) take place. The original conjecture 
of Pettini and coworkers asserts that the two energies 
" correspond^ (let's call this the strong topological hy- 
pothesis). To our knowledge there is only one system 
within the hypotheses of the theorem, the two dimen- 
sional tp 4 model 3], for which the equivalence has been 
numerically established. In other two systems with long 
range interactions, thus out of the theorem hypotheses 
(the mean field XY model and the mean-field k- 
trigonometric model H,@) the equivalence has been ana- 
lytically proved. There are, conversely, analytical results 
for a different model system (mean field <p 4 model 0) 
for which the correspondence does not hold: the energy 
v c at which phase transition takes place is higher than 
the energy vg of the topological singularity (we stress 
here that also in this case the hypotheses of the theorem 
in Ref. || are not fulfilled). 

The latter discussion is of particular importance, not 
only in the context of phase transitions, but also in anal- 
ogy with glassy systems. These systems are characterized 
-at the mean field level- by a dynamical transition tak- 
ing place at a given temperature T MCT (or equivalently 
at energy v MCT ) predicted by mode-coupling theory, and 
a (hypothesized) true phase transition at a lower temper- 
ature Tk (Kauzmann temperature) or energy vk- From 
numerical simulations of Lennard- Jones like systems, one 
observes that the dynamical transition is strictly related 
to the properties of the saddles visited by the system 
[HI Il6| . The concept of "visited saddles" is quantita- 
tively worked out defining a pseudo-potential W—\ W\ 2 
and minimizing it during the dynamic evolution of the 
system, thus obtaining a map M. q : q — > q s associating 
to each equilibrium configuration point q = {qi, . . . , q^} 
a minimum q s of W . When averaged over the dynamic 
trajectory one obtains an energy map: M. : v — ► v s . 
Absolute minima of W (having W(q s )—Q) correspond to 
stationary points (saddles and minima) of V. We note 
that the presence of local minima of W, with W(q s ) ^ 
but small (corresponding to inflection directions in V 
profile) , does not affect the result : the order of vis- 
ited saddles (number of negative eigenvalues of the Hes- 
sian matrix of V) extrapolates to zero at T MCT1 and the 
energy of saddles stays always below the instantaneous 
energy. Moreover the true thermodynamic transition is 
achieved when the number of visited stationary points 
of order zero (minima of V) grows less than exponen- 
tially with the system size (in the glassy terminology 
when the " configurational entropy" or 11 complexity" goes 
to zero). Solvable mean field spin-glass models (p-spins), 



which manifest the same phenomenology of structural 
glasses, corroborate these findings in an analytical way 
|18|. Then, what emerges from glassy systems, is the 
great importance of the underlying stationary points in 
the description of the various transitions (dynamical and 
thermodynamical) taking place in these systems. It is 
worth to note, however, that the consistency of this pic- 
ture beyond mean field is still matter of debate [19j , and 
that the definition of the map M. is not unique also at 
the mean field level, different definitions giving similar 
but not quantitatively equal results |icj|. 

One can argue, in line with the "topological approach 
to phase transitions, that also for non-glassy systems the 
concept of underlying stationary points continues to be 
useful. It is important to emphasize that, in the study 
of the glass transition, is the discontinuity of the av- 
erage density number of underlying stationary points 
that marks the dynamical transition at T MCT . Driven 
by this observation, we recently proposed that the map 
M. : v — > v s has to be applied in order to spot the phase 
transition, i.e., if a topological discontinuity exists at en- 
ergy vg , the phase transition is expected at an energy v c 
such that M(v c )—vg. This has been proved to work (at 
least approximately) in those cases (e.g. the mean field 
<f 4 model) where the original "strong topological hypoth- 
esis" (i.e. coincidence between v c and vg) failed. It is 
worth to point out that those cases where it has been 
proved that v c = vg do not constitute counterexamples 
for the application of the map A4, as in all these cases 
it turns out that v c is a fixed point for the map, i. e. 
A4(v c )=v c . In conclusion, for all the cases investigated 
so far 

0I1HG3, 

it results that whenever a phase tran- 
sition (including also "dynamic" transitions as the glas s 
transition in LJ liquids 0,^1 an d p-spin systems |1S|1 
is present at a certain energy v c , this transition is sig- 
naled by a discontinuity in the topology, specifically in 
the Euler characteristic or in the complexity, at an energy 
vg such that M(v c )=vg. At variance with the original 
(strong) topological hypothesis, where it was supposed 
that v c = vg, we will refer to the latter conjecture as 
weak topological hypothesis. Note that the weak topologi- 
cal hypothesis, at variance with the strong topological hy- 
pothesis, depends on the statistical measure, as the map 
M. : v — > v s is defined through an average over the dy- 
namical trajectory (or, equivalently, over the statistical 
measure). We will discuss this point in detail in the fol- 
lowing. 

Two recent papers 0, 0] addressed the question con- 
cerning the relationship between phase transitions and 
topology in one dimensional models. Kastner stud- 
ied two versions of a solid-on-solid model, one showing a 
phase transition at finite temperature and the other not; 
he found that both models exhibit the same topologi- 
cal change, thus concluded towards an "unattainability 
of a purely topological criterion for the existence of a 
phase transition" Grinza and Mossa ^3 consid- 

ered the Peyrard-Bishop model [2(], 0] , which exhibits 
both a phase transition and a change in the topology, 
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but in this case v c and vg are not coincident [22|. These 
papers contributed to extend the analyzed cases for the 
understanding of necessary and sufficient conditions for 
the topological hypothesis. However, they seem to reach 
contradictory results, one supporting and the other falsi- 
fying the topological hypothesis, even if the investigated 
models share many similarities. 

The aim of this work is to try to clarify this appar- 
ent inconsistency with regard to what has been discussed 
above. In particular, we reanalyze the model investigated 
by Kastner and by Grinza and Mossa. As a result of 
this study i) analyzing the Peyrard-Bishop Model (PB) 
[20I I2H ] we numerically show that the two energies, al- 
though different v c ^ vg, satisfy the 11 weak topological 
hypothesis" , i.e. A4(v c )=vg. We also investigate a slight 
modification of the PB model (allowing non-confined mo- 
tion of the variables), where we are able to study the same 
quantities in absence of a thermodynamic phase transi- 
tion at finite temperature, again finding results in agree- 
ment with the "weak topological hypothesis" . ii) Analyz- 
ing the Burkhardt model we introduce a further parame- 
ter defining the position of the pinning potential, that can 
be moved from the origin, i.e. fully confining potential, to 
infinity, fully non-confining potential. We found that the 
phase transition actually exists for all the position of the 
pinning potential, and its critical temperature goes con- 
tinuously to infinity as the pinning potential position goes 
to infinity. Moreover, we found that also the generalized 
Burkhardt model falls into the class of systems that sat- 
isfy the " 'weak topological hypothesis" , i.e. M.{v c )—vg. As 
the position of the pinning potential is moved toward in- 
finity, the energy v s (T) of the "underlying saddles" tends 
to reach the energy vg at higher temperature; when the 
nonconfining limit is reached, v s (T) < vg for all T. The 
topological singularity is visited only for T — > 00 and this 
is the reason why the phase transition is not observed, 
i.e. T c — > 00. 

The present findings support the idea that also in the 
case of one dimensional models, the relevant topological 
quantity related to a phase transition is obtained from 
underlying stationary points obtained from the map M. . 
As we already noted, the choice of the map M. is not uni- 
vocal: we have chosen the one obtained through W (with 
some ad hoc modifications) , however different choices are 
possible (we mention here, for example, the map obtained 
using Euclidean distances in configuration space |iclll8j h 
The robustness of this conclusion with respect to the pos- 
sible different choices of the map M. is still an open ques- 
tion which goes beyond the scope of this paper. 



I. THE MODELS 

The one dimensional models we study are all defined 
by the Hamiltonian H = J^iLiPi + V({q}i=i...N) 
(m is the mass of each particle), where V is the poten- 
tial energy. We consider two different classes of models. 
The first one, introduced by Burkhardt [2j| as a model 
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FIG. 1: Sketch of the on-site pinning potential Vp 1 \q) for a 
given choice of the control parameter L/R. 



for localization-dclocalization transition of interfaces, is 
defined by the potential energy V^: 



JY 



N 



V^({q}i= 1 ... N )=J2K\ qi+1 -q i \+Y / V p ^(q i ) , (I) 



where K measures the strength of the force between 
neighboring pairs, V p (1 \q) is the on-site pinning po- 
tential, and periodic boundary conditions are assumed 
<7tv+i = 31 • We chose for V p (q) the following form 



K (1) (<z) = 




for q < 
for < q < L 
for L < q < L + R 
for q> L + R 



(2) 



that generalizes the original form in Ref. [23 introduc- 
ing a parameter L that gives the position of the pinning 
potential (a square well of depth C/ and width R, see 
Figure ^| from the edge of the system. The case with 
L = coincides with the original Burkhardt confining 
model, while the non-confining case is retrieved in the 
L — > 00 limit. 

The models of the second class are defined by the the 
potential energy and of the form: 



N 

V^({qh =1 ... N ) = £ y (q i+1 qi f + £ v^\ qi ) . 

i=l i=l 

(3) 

We consider two different versions of this model, one 
defined by the on-site Morse potential, introduced by 
Peyrard and Bishop as a simple model for DNA thermal 
denaturation [2(], (PB model) 

V^(q) = U {(e-<i/ R -lf-l}; (4) 

the other is a symmetric version of the former (SPB 
model) 



v m {q) = Uo{{e -\i\/n- 1 f-i ] 



(5) 
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FIG. 2: Plots of the on-site pinning potentials V p (q) (full 
line) and Vp 3 \q) (dashed line). 



a slight modification of PB model that allows for a non- 
confined motion of the variables (see Fig. EJ). We note 
that the introduction of the modulus in the Eq. [5] does 
not introduce discontinuities up to the second derivative 
of the potential. The quantities U and R determine re- 
spectively the energy and the length scales of the on-site 
potential (in the following all quantities will be reported 
in U and R units). We further chose m=l. In all the 
three cases a parameter of the Hamiltonian is related to 
the strength of the inter-particles interactions (K) and, 
for the case of a second parameter is the position of 
the pinning potential L. 

Specifically, the relevant quantity defining the relative 
weight of the on-site with respect to interparticles poten- 
tials is the dimensionless ratio £=KR/U or £_—KR 2 /2U 
for the potential models (1) or (2,3) respectively, while 
the position of the pinning potential is given by £ = L/R 
for the potential model (1). 

The generalized Burkhardt model has been treated an- 
alytically, while the Peyrard-Bishop models arc studied 
numerically. In the latter cases we performed isothermal 
molecular dynamics simulations using Nose- Hoover ther- 
mostat at different temperatures for systems with iV=500 
degrees of freedom with periodic conditions qjv+i=9i- 
We studied different values of the control parameter £ , as 
an example the results are reported for £=0.05 and 0.5, 
all the other £- values give results in qualitative agreement 
with the two reported examples. 



II. BURKHARDT MODEL 

A. Thermodynamics 

The thermodynamics of the Burkhardt model is known 
since many years [23j for both the £ = and £ = oo cases. 
The method of solution is briefly outlined below. 



The determination of the thermodynamic of systems 
described by a potential function of the form 



JY 



N 



V({q} i= i... N ) = ^,K\q i+1 - qi \ + ^V p ( qi ) , (6) 



i=l 



i=l 



i. e. similar to the case (1) (Eq. QJ, goes through the 
exploitation of the transfer matrix technique. Indeed, 
the configurational partition function Z is given by: 



Z N = I dq\... dq N e 



-/3V(M,=i...iv) 

that, defining the transfer "matrix" 

T(x,y) = e -P K \ x -v\ e -W P (x)+v P (y)]/2^ 

can be written as: 

N 

Z N = / dqi...dq N JJ T( qi , q i+1 ) , 

i—l 



(7) 



(8) 



(9) 



recalling that qN+i = q\- With this notation, the (con- 
figurational) free energy density 



in the thermodynamic limit is promptly written as 
1 







log(max {A}) 



(10) 



(11) 



where A is the set of eigenvalues of the transfer matrix, 
i. e. the eigenvalues of the integral equation 



dy T(x,y)<f>(y) = \<f>(x). 
The latter equation, with the substitution 

turns out to be 



(12) 



(13) 



dye-P^*-^ e~ f)V ^il){y) = \il){x). (14) 



The next step is performed by noticing that the operator 
[—d 2 /dx 2 + (3 2 K 2 ] applied to exp(— [3K\x — y\) produces 
a delta-function: 



o -0K\x-y\ _ 



20K6(x-y), (15) 



thus by applying the previous operator to the integral 
eauation ll4l it can be transformed in a Schroedinger like 
differential equation: 



2(3K 



-PV p (x) 



ip(x) = 0. (16) 
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This equation must be solved with the conditions that i) 
the "eigenf unction" ip(x) was normalizable, and, ii) the 
boundary condition (implicit in EclITUi ip'(0)/ip(Q) = (3K 
was fulfilled. In summary, the calculation of the thermo- 
dynamic of system defined by the potential energy of the 
form in Eq.^is reduced to the solution of a Schroedinger- 
like differential equation and, in particular, to the finding 
the largest eigenvalue of the original integral eauation ll4l 
In general, as the eigenvalues are continuous and smooth 
function of the parameters (among which the tempera- 
ture), no phase transitions are expected unless the two 
largest among them cross each other. 

B. The C = case. 
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Let us now apply the procedure to the potential func- 
tion in Eq. for the case C — 0. We do not report the 
details of the calculation, as they are based on standard 
techniques for solving Schroedinger equation in Quan- 
tum Mechanics [24j ; in summary the " eigenvalues" A are 
determined by the equation: 



2(A) = (3K 



(17) 



with 



m _ fi(P, Q) sin (QR) - UP, Q)) cos (QR) 
Z[ ' fa(P, Q)) sin (QR) + fi(P, Q) cos (QR) ' 1 ' 

having defined 

h(P,Q) = Q 2 (19) 
f 2 (P,Q) = PQ 

h(P,Q) = P 
f 4 (P,Q) = Q 

and 



(20) 



P(X) = \ I3 2 K 2 - 



2(3K 
A ' 



The only possibility for the function z(X) to be real 
(condition required for Eq.Eh have solution) is that P 
was real (if Q become imaginary z(X) is still real), thus 
it exists a solution to Eq. El only if P is real. There- 
fore, when P vanishes, the eigenvalues A disappear (more 
specifically, disappear the eigenvalues of the discrete 
spectrum, and only those of the continuum spectrum re- 
main) , and the (configurational) free energy is discontin- 
uous. For each temperature, the condition P(X) = is 
fulfilled for a "critical" A, given by: 



(21) 



Thus, the equation for the largest eigenvalue at the "crit- 
ical" point is given by z(X c ) — (3 C K, or 



y / e^ u " - 1 tan p c KRy/eP°V' 



1 = 1, 



(22) 



FIG. 3: Full line: critical temperature (in reduced units 
kBT c /Uo) as a function of the control parameter £ from Eg. 1231 
for the £=0 case. Dashed line: temperature Tj (in reduced 
units) at which the "underlying saddle" jumps from minimum 
to saddle as a function of £ (see Sec. V B). 



which gives us the required equation for the critical (in- 
verse) temperature (3 C . This equation can be rearranged, 
introducing the control parameter £ = KR/Uq, as: 



1 



1 



PcUo Ve^o _ i 



: arctan 



1 



y/eP° u ° - 1 



(23) 



The plot of the critical temperature (in reduced units 
ksT/Uo) as a function of £ is reported as full line in 
Fig. 3. 

The phase transition is of the localization- 
delocalization type. The particles, kept together 
by the K\x — y\ term of the potential, for T < T c are 
pinned close to the square well, while, for T > T c are 
delocalized in the g-axis. 

A simple calculation leads to the value of the critical 
energy v c (the equilibrium energy v(T) at the transition 
point v c = v(T c )). From Eq. 1111 we have 



v(T) 



0(0/) A'(/3) 



d(3 



A(/3) 



(24) 



where X(f3) is the solution of Eq.EI Close to the critical 
point, AOS) = 2/f3K, thus X'(f3)/X(f3) = 1/(3 and 

v c = k B T c (25) 

independently from the value of £. 

C. The C case. 

The calculation for the case of generic £ values is quite 
similar to the previous one. Also in this case, the eigen- 
values A are determined by an equation like Eq. 1171 
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FIG. 4: Critical temperature (in reduced units ksTc/Uo) as 
a function of the control parameter £ from Eq. 1281 for the 
indicated £ value: £ = (full line), 0.5 (dashed line), 1.5 
(dotted line), 3.5 (dash-dotted line) and 7.5 (dash-dot-dotted 
line). 

z(X) — (3K, with z(A) again given by Eg. 1181 and with 
the f n {P, Q) functions (n—1, . . . , 4) given by 

fi{P, Q) = P [Q 2 cosh(PRC) - P 2 sinh(Pi?C)] 

h{P,Q) = P 2 Qexp(PR() 

fa(P,Q) = [P 2 cosh(Fi?C) - Q 2 sinh(Pi?C)] 

U(P,Q) = PQexp(PRC). (26) 

Obviously, Eas. 1261 recover Eqs. |201in the £ — > limit. 
The same considerations on the reality of P(A) reported 
above apply to Eq. Thus the condition P(A) = 
define the critical value of the eigenvalue, A c = 2/ j3K, and 
the equation for the critical temperature (z(X c ) = (3 C K) 
becomes: 

y/ e 0cUo _ i S i n (j3 c KR^Je^ u ° - 1 ) x (27) 
x [cos (j3 c KR^Je^ u ° - 1 ) - 

-/3 c KRC\Ze^ u " - 1 sin (p c KRyfeW° - 1 = 1. 

Similar to the £ = case, this equation can be rear- 
ranged, introducing the control parameter £, as: 

_L 1 / 1 1 \ 

^ " Mo y/eW- - 1 arCtan We^-1 T+^oCC J ' 

At variance with Eq. 1231 this equation cannot be cast in 
the form £ = £(/3 c ), thus it must be solved numerically to 
plot the critical temperature as a function of the control 
parameter £. This plot is reported in Fig. 0] for different 
values of (. 

As can be observed in Fig. 0] on increasing £, i. e. on 
displacing the position of the square well towards high 
value of the coordinate, the critical temperature, for a 



FIG. 5: Critical temperature (in reduced units kBT c /U ) as 
a function of the control parameter £ from Eq. 1281 for the 
indicated £ value: £ = 2 (full line), 1 (dashed line), 0.5 (dotted 
line), 0.25 (dash-dotted line). 

given £ value, increases, expanding the amplitude of the 
"cold" (localized or pinned) phase. This can be better 
seen in Fig. [5] where the ( dependence of the critical tem- 
perature is reported for some values of £. We conclude 
this section noticing that the phase transition actually ex- 
ists for all the value of £, and in the limit of £ — > oo, the 
critical temperature goes without discontinuities to infin- 
ity. Therefore, we are lead to conclude that the model 
investigated in Ref. 01 to demonstrate the unattainabil- 
ity of a purely topological criterion for the existence of 
a phase transition is a "borderline" model, in which the 
phase transition can be thought to be present at "T infin- 
ity" (even though the precise meaning of this statement 
is not well defined). Then, to the same topology (as we 
will see in the next section, the topology of the poten- 
tial function in Eq. [2 does not depend on the value of C) 
always corresponds a phase transition. 

To discuss the question of the coincidence (or not) of 
the critical energy with the topological discontinuity, we 
need to calculate v c (£,0- Following the same argument 
reported for the case £ = we conclude that the critical 
potential energy depends on £ and £ only through T c : 
v c = ksT c . As an example, in Fig.[S]we report the caloric 
curve v(T) as a function of the inverse temperature for 
different C value and for £ = 1. For all the C values, on 
the low-/? side the curves end at the points [fi c , v c ); these 
points are aligned along the v c (j3 c ) line (thick dotted line) 
given by v c (fi c ) = 1/(3 C . 

D. Topology 

The analysis of the topological properties of the 
Burkhardt model is reported by Kastner in Ref. [Tl| . 
He analyzed only the two limiting cases of confining and 
non-confining models, corresponding in our notation to 
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FIG. 6: Inverse temperature dependence (in reduced units 
PU ) of the equilibrium potential energy for £ = 1 and for 
the indicated £ value: £ = (full line), 0.5 (dashed line), 1.5 
(dotted line), 3.5 (dash-dotted line) and 7.5 (dash-dot-dotted 
line). The thick-dotted line represents the /3 dependence of 
the potential energy in the high temperature phase. 



£=0 and £=oo respectively. He found that a topology 
change is present in both cases, even if not really equal in 
"strength". The value of the potential energy at which 
the topological change appears is vg=0, irrespective of 
the considered model. One can easy generalize the above 
analysis to the general case with arbitrary £, and con- 
clude that the topological change is always located at 
energy vg=0. It is worth noting that the energy at which 
topological change appears is lower than the thermody- 
namic transition energy: v c > vg (see Fig. |SJ). We will 
further discuss this issue in Sec. V, after having described 
the thermodynamics and topology of the PB and SPB 
models. 



III. PEYRARD-BISHOP MODEL 

A. Thermodynamics 

The thermodynamics of the Peyrad-Bishop model (de- 
fined in Eq.s |3 can be studied using transfer matrix 
techniques, as the Burkhardt model described in the pre- 
vious section. However, in this case approximated meth- 
ods have to be considered in order to obtain a corre- 
sponding Schroedinger like differential equation. In the 
region £ 3> 1 and temperature window U Q *C k B T <C £JJ D 
the classical statistical mechanics problem is mapped to 
the quantum Morse oscillator problem [2(|. Simi- 
larly to the case of Burkhardt potential, the presence of 
a second order phase transition for the Peyrad-Bishop 
model is signaled by the bounded-unbounded transition 
of the lower state in the corresponding quantum prob- 
lem. In the above range of £ and T, Peyrard and 
Bishop obtained an analytical expression for the tran- 




FIG. 7: Temperature dependence (in reduced units k B T /U ) 
of the equilibrium potential energy v (full symbols) and en- 
ergy v s of underlying saddles (open symbols) for the PB 
model defined by Eq. 0]with £=0.05 (upper panel) and £=0.5 
(lower panel). Also indicated in the hgure are the values of 
topological change energy vg (horizontal dotted lines), phase 
transition energy v c (horizontal dot-dashed lines) and transi- 
tion temperature T c (vertical full line) for £=0.05 (yg /U o =0, 
Vc/Uo ~ 0.61 and k B T c /U ~ 1.22) and for £=0.5 (v e /U o =0, 
Vc/Uo — 1-59 and kBT c /U — 3.20). Dashed lines are the 
T-dependence of the potential energy in the high T phase: 
v(T)=k B T/2. 



sition temperature fcsX' c /[/ =4v / £ and transition energy 
v c /U =kBT c /2U =2y/^. For generic (£,T) values, only 
numerical results can be used to infer the existence and 
location of a phase transition. In Fig. 0we report the 
temperature dependence of the potential energy per par- 
ticle v — V/N (full symbols) of the PB model for two 
different values of £: 0.05 (upper panel) and 0.5 (lower 
panel) . Also reported in the figure are the energy v c (dot- 
dashed line) and temperature T c (full line) of the phase 
transition point: v c /U — 0.61 and kBT c /U — 1-22 for 
£=0.05, Vc/Uo ~ 1.59 and k B T c /U a ~ 3.20 for £=0.5. 
Dashed lines are the T-dependence of the potential en- 
ergy in the high T phase: v(T)=k B T/2. 



B. Topology 

The topology of the Peyrard-Bishop model is studied 
in the paper of Grinza and Mossa [T3I - A topological 
change is found at the energy value vg=0, corresponding 
to a topological change in the hypersurfaces T, v varying 
v: from a close hypersurface for v < vg to an open one 
for v > vg 12]. In Fig. 0the value of vg is indicated by 
an horizontal dotted line. We note that, also in this case, 
the topological discontinuity is lower in energy than the 
thermodynamic one: v c > Vg. 
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FIG. 8: Temperature dependence (in reduced units fcsT/!7 ) 
of the equilibrium potential energy v (full symbols) and en- 
ergy v s of underlying saddles (open symbols) for the SPB 
model defined by Eq. 0with £=0.05 (upper panel) and £=0.5 
(lower panel). Also indicated in the figure is the value of 
the topological change energy ve/U o —0 (dotted line) for both 
cases. 



IV. SYMMETRIC PEYRARD-BISHOP MODEL 

A. Thermodynamics 

The Symmetric Peyrard-Bishop model defined by Eq.s 
13 [S] does not exhibit phase transition at finite T. This 
can be view from the fact that there is always a bound 
state in the corresponding quantum problem, in analogy 
with the non-confined Burkhardt model |23l |27| . In Fig. 
|H1 the same quantities as in the PB case are reported 
for the SPB model: energy v (full symbols) for £=0.05 
(upper panel) and £=0.5 (lower panel). It is evident in 
this case the absence of a phase transition at finite T (in 
the T-range investigated) . 



B. Topology 

Following a similar argument as in Ref. |Tlj . one can 
see that also in the SPB case one has a topological change 
at exactly the same energy level as in the PB model vg=0 
(even if not identical in stren gth to the previous one). We 
refer to the papers in Ref.s [HQ2 for 

a more detailed 

discussion of the topology. In Fig. |H1 the value of vg is 
indicated by an horizontal dotted line. 



V. UNDERLYING SADDLES 

In this section we study the properties of the stationary 
points visited by the systems. The concept of "underly- 



ing saddles" was first introduced in the study of glassy 
disordered systems 0, 0, to better understand the 
topological counterpart of the dynamic transition taking 
place in these systems. Recently, it has been applied 
also in the analysis of models that exhibit thermody- 
namic phase transitions, in order to emphasize the role 
of topological changes at the "underlying saddles" energy 
in driving the phase transition |8j, |9|, llO( . 

Here we apply the same methodology to investigate the 
one dimensional systems introduced before. Let start 
with the models having a continuous potential energy 
function, the PB and SPB models, which allow for the 
usual definition of stationary points. At the end of the 
section we will extend the argument to the discontinuous 
case of Burkhardt model. 



A. Peyrard-Bishop and Symmetric Peyrard-Bishop 
models 

There are only two stationary points in the potential 
energy hypersurface of both models: a minimum located 
at qi=c[2=. . . =(/Ar=0 and a saddle (with degenerate Hes- 
sian matrix) at qi=q2=- ■ ■ =<zat=oo 01 • ^ n order to as- 
sociate one of the two stationary points to each instan- 
taneous configuration of the system, we used a similar 
trick as in the anal ysis of glassy systems or mean- 
field models 0, |j| TIo| . In the latter one minimized 
the pseudo-potential PF=|Vy| 2 during the dynamic evo- 
lution at different temperatures, so introducing a map 
from equilibrium energy levels to saddles energy levels: 
M. : v — ► A4(v) = v s . Due to the peculiarity of the 
present models, where the saddle point is "infinitely" far 
from each equilibrium configuration, we decided to apply 
the W minimization method in a two steps procedure: i) 
first we minimized the Wi n t quantity defined using the 
interaction potential part of V, Wi n t=| Wi„ t | 2 , where 
Knt=Ei=i t(*+i " *) 2 ! **) tncn wc minimized the W p 
defined using the on site potential W p =| Wp 2 ' 3 ' | 2 . This 
procedure ensures that the point reached is a true sta- 
tionary point, i.e. the minimum or the saddle. Obviously, 
this is a quite arbitrary definition of basins of attraction 
of stationary points. As said in the introduction, the ro- 
bustness of the results with respect to the possible choices 
of definition of a saddle basin of attraction is still an open 
problem. 

In Fig. 13 the temperature dependence of the energy 
v s (open symbols) of underlying saddles is shown for the 
case £=0.05 (upper panel) and £=0.5 (lower panel) in the 
PB model. The remarkable fact is that at T c (vertical full 
line in Fig. [JJ the identity v s =vg holds. The map M. (v) is 
shown for PB model (open symbols) in Fig. Elfor the two 
cases £=0.05 (upper panel) and £=0.5 (lower panel). One 
observe that, as before pointed out, one has M(v c )—vg 
for both £ values. The fact that v s (T) in Fig. as well 
as M.(v) in Fig.^l has a "smooth" transition between its 
low T (or v) and high T (high v) regions is most likely 
due to a finite size effect (iV=500 here) and both v s (T) 
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FIG. 9: Map M : v — » v s defined minimizing the pseu- 
dopotential W=\W\ 2 in the PB and SPB models for £=0.05 
(upper panel) and £=0.5 (lower panel). Also reported are 
the corresponding vg (dotted lines) and v c (dot-dashed lines) 
for £=0.05 (i-9=0, v c ~ 0.61) and £=0.5 (v g =0, v c ~ 1.59), 
evidencing the identity A4(v c )—vg. 



and M(v) will probably tend towards a step function in 
the thermodynamic limit. The previous finding indicates 
that the relevant quantity to consider when we are look- 
ing for topological changes related to a phase transition is 
the underlying stationary point energy, obtained trough 
a map from the critical level v c . It is worth noting that 
the map M. is constant (A4(v)—vg) for a broad range of 
values, also below v c , at variance with other cases where 
around the transition point the properties of visited sad- 
dles change 0, El 0- One can conjecture that the 
flatness of M.(v) is a pathology of these one-dimensional 
models, that have a number of stationary points that is 
not extensive in N (actually there are only 2 stationary 
points). 

In Fig |H1 we report the same quantities v s as before 
(open symbols), now for the SPB model, with £=0.05 
(upper panel) and £=0.5 (lower panel). In this case no 
phase transition is present, and indeed the topological 
singularity is never visited, v s (T) < vg for each finite 
temperature (T < oo). 



B. Burkhardt model 

To apply the analysis of the previous section also to 
the Burkhardt model, one has to find a suitable defini- 
tion of "saddles" and of "basin of attraction of a sad- 
dle" for a discontinuous potential. One possibility is 
the following: we first minimize the interaction poten- 
tial Vi n t = Yli=i K\<li+i ~ which is equivalent to 
put all the qi equal to the center of mass coordinate 
q = N^ 1 J^iQi- ^ Q hes m the well of the potential, 



i.e. q E [L.L + R], we will associate the "minimum" to 
the initial configuration, otherwise we will associate it to 
the "saddle" (we use this terminology by analogy with 
the PB model) . It is clear that the average energy of the 
"underlying saddles" is simply the average of the on-site 
energy of the center of mass coordinate, 



T ■ 



(29) 



In the thermodynamic limit the center of mass q is 
peaked around its mean value and then we can substi- 
tute the right hand side of Eq. l2~9l with V p ((q)), a 
quantity that can be explicitly computed. To determine 
(q) we can use the distribution probability |0(x)| 2 , where 
(j)(x)^e~ /3Vp< ' x ^ 2 ^(x) and ip(x) is the eigenfunction of the 
transfer matrix operator corresponding to the maximum 
eigenvalue [23| (see Sec. II A). We note that the sad- 
dle energy v s (T) is a step function, equals to the mini- 
mum energy — U when (q) lies inside the square well and 
equals to the saddle energy otherwise. The tempera- 
ture Tj at which the visited "underlying saddle" jumps 
from minimum to saddle is shown in Fig. [3] (dashed line) 
as a function of the parameter £ for the £=0 case. It is 
worth noting that the temperature Tj lies always below 
the thermodynamic transition temperature T c (in anal- 
ogy with the PB model, see Fig. 0). The same happens 
for all values of £. Therefore, also for the Burkhardt 
case, at the transition temperature T c the "underlying 
saddles" lie at an energy equal to the topological discon- 
tinuity energy vg, i.e. M{v c )=vg. 



VI. CONCLUSIONS 

Studying two particular one dimensional models dis- 
cussed in the recent literature 

[HG3 (Burkhardt model 
in the confining and non-confining version, Peyrard- 
Bishop model and its non-confining counterpart), we 
have focused on the relationship between phase transi- 
tions and topological changes, recently proposed in the 
literature @,H,SI3- I n these models, a topological sin- 
gularity at a given energy value vg(=0) is always found; 
however, i) in the confining version a phase transition is 
found but the critical energy is v c > vg ii) in their 
non-confining version there is no phase transition at any 
finite temperature [Tlj |. 

These results generated confusion as i) was interpreted 
as a confirmation of the strong topological hypothesis of 
Pettini et al. while ii) was considered as an evidence 
for the unattainability of a purely topological criterion for 
detecting phase transitions, although demonstrated only 
for the particular non-confining one dimensional models. 

Exploiting the concept of 11 underlying stationary 
points" defined through a generalization of the meth- 
ods used in the glassy literature (minimization of the 
pseudopotential W / =|VF| 2 ), we have defined a map M : 
v — > v a from energy level v of V to stationary points, 
with energy v s . We have shown that: i) in the con- 
fining case, where the phase transition is present, one 
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has M.(v c ) = vg, in agreement with the weak topologi- 
cal hypothesis; ii) in the non-confining case, where the 
phase transition is not present at finite temperature (as 
the transition temperature goes continuously to infinity 
when the confining wall is removed) the energy of the 
underlying saddles is always below the topological singu- 
larity, i.e. v s (T) < vg , VT; the singular point vg is indeed 
visited for T — > oo, consistently with the observation that 
the critical temperature is "infinite" in the non-confining 
case. 

The weak topological hypothesis appears as a possible 
framework to fit the results that recently appeared in 
the literature on all the different models investigated so 
far. Within this hypothesis three different scenarios are 
possible: 

1. If there is no topological singularity vg, a phase 
transition is not possible; this is consistent with the 
hypothesis of Pettini et ai: topological singularities 
are necessary conditions for a phase transition to 
take place. 

2. If there is a topological singularity at energy v$, a 
phase transition is also present if and only if there 
exist a temperature T c such that v s (T c ) = vg (or 
equivalently an energy v c such that M(v c ) — vg). 



The above findings seem to indicate that, at least for the 
particular models investigated, a sufficiency criterion for 
the phase transition to take place requires the introduc- 
tion of a statistical measure: thus, we believe that the 
statement of Kastner [Tlj concerning the unattainability 
of a purely topological criterion for detecting phase tran- 
sitions is indeed correct, even though in Ref. |nj it has 
been derived using a "borderline" model (see Section II 
C). 

Let us conclude with two remarks: i) as already stated, 
the definition of the map M. is not unique, different defi- 
nitions giving (slightly) different results. Thus, the weak 
topological hypothesis contains in its formulation an am- 
biguity and must be regarded only as a practical tool, at 
least at this stage of comprehension; ii) nevertheless, we 
hope that this approach can be of interest for the numer- 
ical investigation of systems of "mesoscopic" size (e.g. 
proteins and large molecules), i.e. such that the num- 
ber of degrees of freedom is not large enough to allow to 
detect the presence of a phase transition using standard 
techniques. 

We thank M. Pettini and M. Kastner for helpful com- 
ments and suggestions. 
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